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Abstract 

The behavior of rj and rj' in hot strange quark matter in weak equihbrium with temperature, 
is investigated within the SU(3) Nambu-Jona-Lasinio [NJL] model. Possible manifestations of 
restoration of symmetries, by temperature or density, in the behavior of r/ and r/' are discussed. 
The role played by the combined effect of temperature and density in the nature of the phase 
transition and meson behavior is also analyzed. 
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1 Introduction 



Understanding the QCD phase diagram is one of the great challenges in the physics of strong 
interactions. It is expected that under extreme conditions (high density and or temperature) phase 
transitions occur leading to a plasma of deconfined quarks and gluons [QGP] [H El E] • While the 
phase transition at zero chemical potential and finite temperature is accepted to be second order 
or crossover, there are indications that the phase transition with finite chemical potential and zero 
temperature is first order. Experimental and theoretical efforts have been done in order to explore 
the /i — T phase boundary. Recent Lattice results indicate a critical "endpoint", connecting the 
first order phase transition with the crossover region at = 160 ± 35MeV , he = 725 ± 35MeV 
|3]. Besides restoration of chiral symmetry and deconfinement, the rich content of the QCD phase 
diagram has been recently explored in other directions. For instance, in high density matter, for 
which neutron stars provide an excellent laboratory, a color locked phase with kaon condensates is 
expected to occur. 

Much attention has also been paid to the question of which symmetries are restored. In the 
limit of vanishing quark masses the QCD Lagrangian has 8 Goldstone bosons, associated with the 
dynamical breaking of chiral symmetry. In order to give a finite mass to the mesons chiral symmetry 
is broken ab initio by giving current masses to the quarks, therefore in the high density/temperature 
regime an indication of the restoration of chiral symmetry would be that the mesons became again 
almost Goldstone bosons. However, the mass of r]' has a different origin than the masses of the other 
pseudoscalar mesons and it cannot be regarded as the remanent of a Goldstone boson. The problem 
of the non-existence of the ninth Goldstone boson, predicted by the quark model, was solved by 
assuming that the QCD Lagragian has a Ua{^) anomaly. The mass of the rj' is obtained by explicitly 
breaking the Ua{^) symmetry, for instance by instantons. 

The study of the behavior of mesons in hot and dense matter is an important issue since they might 
provide a signature of the phase transition and give indications about which symmetries are restored. 
There is a long standing question ^ about whether only SU(3) chiral symmetry is restored or both 
SU(3) and Ua{1) symmetries are restored. The behavior of i]' in medium or of related observables 
like the topological susceptibility might help to decide between these scenarios. A decrease of the rj' 
mass in medium could manifest in the increase of the r]' production cross section, as compared to 
that for pp collisions |i6j . 

This paper is devoted to investigate the in medium behavior of the rj and r]' mesons in hot and 
dense matter in the framework of the SU(3) Nambu-Jona-Lasinio [NJL] [7^ model including the 't 
Hooft interaction which breaks the Ua{^) symmetry. At variance with the other pseudoscalar SU(3) 
mesons there is not much information about the behavior of these mesons in medium. Concerning 
the 7] meson, the recent discovery of mesic atoms might provide useful information in this concern 

E 

The NJL model is an effective quark model where the gluonic degrees of freedom are supposed 
to be frozen and, besides its simplicity, has the advantage of incorporating important symmetries 
of QCD, namely chiral symmetry. Since the model has no confining mechanism, several drawbacks 
are well known. In what concerns our present problem, it should be noticed that the t]' mass lies 
above the qq threshold and the model describes this meson as qq resonance, which would have the 
unphysical decay in qq pairs. For this reason, it has been proposed that in order to investigate 
the possible restoration of symmetry, instead of m^/, the topological susceptibility x? a more 
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reliable quantity in this model, should be used p. 

The behavior of SU(3) pseudoscalar mesons, in NJL models, with temperature has been studied 
in |10t ITTj . Different studies have been devoted to the behavior of pions and kaons at finite density 
in flavor symmetric jT2| or asymmetric matter El Ej • The combined effect of temperature and 
density on kaons in quark matter simulating symmetric nuclear matter has been investigated in 

In this work we investigate the phase transition and the behavior of rj, rj in strange quark matter 
in weak equilibrium at zero and finite temperature. We present the model and formalism in the 
vacuum (Section 2) and at finite density and temperature (Section 3). In Section 4 we discuss 
present our results for the phase transition in the T — p plane and in Section 5 we present our results 
and conclusions for the r] , rj' behavior. 

2 Model and Formalism 

We will consider in this paper the U{1)a anomaly in the three flavor Nambu - Jona - Lasinio 
model. The Lagrangian we employ here can be used as an effective theory of QCD and has the 
following form 



L = g(,a-7-m)g+f ^[(gA'^g)V(g(z75)A'^gr 

^ a=0 

^ a=0 

+ 5(D[det[g(l - 75)^] + det[g(l - 75)^] . 



(1) 



Here q = {u, d, s) is the quark field with three flavors, Nf = 3, and three colors, = 3. A'^ are the 
Gell - Mann matrices, a=0, 1, ... ,8, A'^ = -y/fl- The explicit symmetry breaking part contains 
the current quark masses rh = diag(m„, m^, m^). Note, that m^ ^i ~ 10 MeV, ~ 100 MeV at the 
scale of QCD (~ 1 GeV). The last term in (^Q) is the lowest dimensional operator and that has the 
SUl{3) ® SUji{3) invariance but breaks the Ua{^) symmetry. This term is a reflection of the axial 
anomaly in QCD. The model is fixed by the coupling constants gs,9v,9D, the cutoff parameter A, 
which regularizes momentum space integrals and current quark masses. In order to fix parameters 
we will use the values of pseudoscalar meson masses, decay constants and the quark condensates at 
zero temperature and baryon number density. 

Using the definition of the determinant we can rewrite the Lagrangian (Q) as 



a=0 



a=0 



+ Y [lDahc{qX''q)iqX\)iqX'q) - Da,,{qX'^q){q{t^,)X'q){q{t^,)X'^q) 



(2) 



with summation over a,b,c = (0, 1, 2, . . . , 8). The structure constants Date are coincident with SU(3) 
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structure constants dabc and D, 



^Oab 



--^^ab and Dqoo 



|. Making the shift 
iq\''q)+ < qX^q > , 



(3) 



where < qX'^q > is the vacuum expectation value, and keeping second order terms of {qX"'q), we get 
for pseudoscalar and scalar interactions (for simplicity, we omit vector and axial - vector parts) 



L = q{td-^-m)q+ ^{{qX'^q)Sab{qX'q) + (g(^75)A"g)Pa6(g(^75) A)} 
in terms of the projectors 

Sab = gsSab + goDabc < qX^q > , 

Pab = 9sSab - goDabc < qX^q > ■ 

The hadronization procedure is made via the functional integral 



dqdq exp {iWeff[q,q]} 



(4) 



(5) 
(6) 



(7) 



with the quark effective action 



^eff[q,q] = / dxL{x 



dx 



q {id ■ 'J — rh) q 



+ l{q\^q)Sab{qX'q) + \{q{i^,)X'^q)P,bmi,)X'q) 
After integration over quark fields in (jZj), we obtain the effective meson action 



(8) 



2 

-iTi In 



(9) 



Here the symbol Tr means the summation over discrete indices and integration over continuous 
variables. The fields a"" and y?" are scalar and pseudoscalar meson nonets, respectively. 
The first variation of the action (jH)) leads to the Dyson - Schwinger equation. 



Mi 



rrii 



'^9s < qiqi > -2fi'D < qjqj >< qtqk > 



(10) 



with k = u,d,s cyclic. Mj are constituent quark masses. The equation (fTUj) is a system of three 
coupled equations and the flavor mixing occurs through the coupling constant ga. It is an effect of 
anomaly contribution. The quark condensates in (jlUj) are determined by 



< qiqi >= -^Tr 
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p-Mi 



-iTi [Si{p)] . 



Here Si{p) is the quark Green function. 

To consider the meson mass spectrum, we expand the effective action (jU)) over meson fields. 
Keeping the pseudoscalar mesons only, we write the second order effective action as 



(2) 



e//L 



--09" 

2^ 



P~,' - H,,(P) 
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(12) 



with the polarization operator Ilab{P), which in the momentum space has the form 

' P. ^ , P. 



Iiab{P)=iN,j ^tr^ [^^(P+^)(A%(^75)5,(P-^)(A%(^75) 



(13) 



where now tr/j is the trace over Dirac matrices. In the result, we have obtained an unnormalized 
inverse meson propagator D^^ = P^i^^ — Ilab{P) or 



D 



1 



ab 



^ab 



^ab{P) 



(14) 



which is the 3x3 matrix in the flavor space. We will neglect tt^ — r] and 7r° — r]' mixing terms because 
they are proportional to < quQu > — < qd(ld >■ We will do the approximation of neglecting also this 
contribution in the medium. In the basis of rj — rj' mesons we may find Pab and Ilab{P) 



Poo Pqs 
Po8 Pss 



^ab 



A 



Pss 

—Po8 



-Po8 

Poo 



(15) 



with the determinant A = Pqo-Pss ~ Pqs- Then the inverse meson propagator ()14|) can be presented 
in the explicit form 





Pss — 


AHoo 


-A 


^ —Pos - 


- AHos 




' A B 




-A 


B C 





~-fo8 ~ Alios 
-Poo ~" Allgg 

' A \ 
C 



O 



(16) 



with the orthogonal transformation matrix O. Let us suppose that the matrix O has the form 



O 



cos 9 sin 9 
— sin 9 cos 9 



(17) 



The value of the angle 9 can be fixed from the condition 

2B 



tan 29 



C-A 



which guarantees us a diagonal form of D^^. From (fTTjl and p8|) we may find A and C_ and write 
Dnh as 



A 
C 



o = \o-^ 

2A 







D-,' 



O 



where 



and 



= {A + C)- ^{C - Af + 452 



D^/ = [A + C) + ^{C - Af + 452 
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(19) 



(20) 
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The masses of the rj and t]' meson can now be determined by the conditions 

D-\M,,O)=0, 



D-/(M^,,0) = 0. 



To calculate the coupling constants, we will express Dab directly in the following form 

Dab = A 



A B \ ^ A f C -B 



Here 



EC D\-BA 



A — -Poo-Pss — -Pqs ' 
D = AC-B^ 

= A(l - tr(Pn) + detPdetn) 



and the effective action (IT2l) now is 



with 

= 5 ( -B A ) • 

It has singularities when D equals zero. Making the pole approximation, we expand D over 

D(P2) = D(P2 = M^) + -^ (P2 - . 

The first term equals zero when P^ = (or P^ = M^,). The second term can be rewritten 
rest frame as 

D(p2) ^ (p2 _ j\^2) 

In the case P^ = we have the following expression for the meson propagator 



( 90 

\dPo 

The matrix 



_ 2M^ 1 / C -B 



allows to parameterize the meson propagator in terms of coupling constants gir) with i = 0,8. Using 
(jHUj) . we obtain 



"2 - ^ (32) 



fl'Orjfl'Sr, = -r^^^-N (33) 



From these coupling constants we may calculate 

9mq = \l^9ov + ^^sr, , (34) 
_ [2 2 

The coupling constants gori', dsri', g-q'qg, drj'ss can be evaluated in the same way by the replacement 
Now we consider the properties of rj and t]' mesons. 

First we solve the gap equation (fTn|) . Then the matrices A,B,C,D depend on the solutions of 
Mi and contain the integrals: 

and 

^^^^^^ = '^'IwT'iP'-Mmp+py-Miy ^^^^ 

which are divergent. To calculate the physical observables we introduce the cut - off parameter A. 
Putting the integrals to (j^^ we may calculate M^j. But, since NJL does not confine the quarks and 
the r]' is above the QuQu and qaQd continuum, the meson propagator has complex poles. We 
assume (see details in ^T] ) that the equation for the poles in the propagator has solutions of the 
form: 

Po = - ^iT, (38) 

where M^/ is the mass and F is the width of the t]' - resonance. 

In order to avoid the complexity of the analytical continuation, we calculate the integrals (jHTjl 
in an approximate way, which means that we neglect F that occurs in the denominator of /2*(Po) 
and consider the small value of the width. Since this last integral is complex, a straightforward 
calculation leads to the replacement in A, B, C, D: 

Plf^iPo) ^ [P^Rer^{Po) + PoTlmP^iPo)] + i [P^\mP^{Po) - PoTReP^{Po)\ |p2=m^,- (39) 
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Therefore from the zeros of the real and imaginary parts of the propagator we get the mass and 
the width of r( . 

For our numerical calculations we use the parameter set: 
m„ = mrf = 5.5 MeV, m, = 140.7 MeV, gst^^ = 3.67, qdA^ = -12.36 and A = 602.3 MeV 
that has been determined by fixing the conditions: 

= 135.0 MeV, Mk = 497.7 MeV, = 92.4 MeV and M^, = 960.8 MeV. 
Here is the pion decay constant jTUl E] . 

Also, we have: 

= 514.8 MeV, e{M^) = -5.8°, g^uu = 2.29, g^ss = -3.71 

and 

M^, = 960.8 MeV, e{M^,) = -43.6°, gr^^uu = 13.4, ^^/g, = -6.72. 

Note that the t]' - meson always lies above the quark - antiquark threshold. Since this is an 
artifact due to the lack of confinement it is a resonant state. Nevertheless, we use M^/ as an input 
parameter, because as mentioned in P|, the chiral susceptibility is related with M^/ by Witten - 
Veneziano mass formula ^3 ^] 

= + - 2Mi (40) 

and takes the value of M^i which is closed to the experimental value. In this formula Nf = 3 
and the chiral susceptibility equals x = (177.2MeV)''. The lattice calculations give the number 
X = (175 ± 5MeV)^ 

For the quark condensates we have < quQu >=< qdQd >= — (241.9MeV)^ and < g^Q's >= 
-(257.7MeV)3. 



3 The finite temperature and density case 

The generalization of the NJL model to the finite temperature and density case can be done if 
we introduce the thermal Green function, which for a system of quarks at given temperature T 
and chemical potential /i, is 

S.{S-x',r-r') = ^Ee-^-^-^U .tX , (41) 

where j3 = 1/T, T is the temperature and the sum is made over the Matsubara frequencies Un = 
(2n + l)vrT, n = 0, ±1, ±2, . . ., so that po — ^ ioOn + /^i with a chemical potential /ij. Instead of 
integration over po "we have now the sum over Matsubara frequencies. Ei = yjp"^ + is the quark 
energy. 

Because in the gap equation (fTUI) instead of Si{p) we put the thermal Green functions (HH), the 
current quark masses Mi will now depend on the temperature and chemical potential 

Mi = nii -2gs « QiQi » -2 go « QjQj »« qtqk » (42) 
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where << QiQi » are the quark condensates at finite T and chemical potential /ij. The condensates 
are expressed in the terms of the integral J|(T, /Xj) which is calculated by the substitution of pijl in 



with the Fermi distribution functions 



{nf - ) 



1 + exp [(3{Ei ± /ij 
The integral l2{Po,T,fii) may be found in the same fashion 



(43) 



(44) 



/2X-P0, T, /ij 



27r2 



V 



2 V « 



'47r\ 



4M2 / ^ P. Po. 
1 TT \nj{ — ) — n, — 



(45) 



Using these integrals we solve equations (j22p and (|23|) at given temperature and chemical potential 
taking into account also the relations 



and 



Poo(^, ^i) = gs- -^9d (<< Qu Qu » + « Qd qd » + « qs qs >>) , 

P8s{T,^i) = gs + ^gD{2 «ququ» +2 «qdqd» - « qsqs»] 



Po8(T,/i.) = P, 



80 



3^2 



9d (<< ququ » + « 



» -2 << qsqs >>) 



(46) 
(47) 
(48) 



where 



^oo{Po,T, iXi) 

n88(Po,7', /ii) 

no8(Po,7', /ii) 



2 [<^M«(-Po! fj,i) + Jdd{Po, T, fii) + Jss(Po, T, /ij)] , 
^ [-^™(^o, 7", /ii) + Jdd{Po, T, fii) + 4Js,(Po, T, /ij)] , 
[Juu{Po, T, fii) + Jdd{Po, T, fii) — 2Jss{Po, T, fii)] 



J,,(Po,T,/i,)=4(Ji + 5L/^0. 



(49) 
(50) 

(51) 
(52) 



4 Phase transition 

To consider the density effects, in accordance with ^3^1^1^112011^; we concentrate here 
on the case of asymmetric quark matter with strange quarks in chemical equilibrium maintained by 
weak interactions and with charge neutrality, by imposing the following constraints on the chemical 
potentials of quarks and electrons and on its densities: 
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Figure 1: Left panel: Pressure T = MeV (1st order phase transition) and T = 56 MeV (2nd order 
phase transition); Right paneh Pressure as function of the density and temperature. At T < 56 MeV 
we have a first order phase transition. After that, the phase transition becomes of second order. 



fJ-d = fJ's = fJ'u + fJ^e and -Pu- -{Pd + Ps) - Pe = 0, (53) 

with 

P^ = ^,i^^^-Mrf'e{p^-M^) and Pe = (54) 

As discussed by several authors, this version of the NJL model exhibits, at T = MeV a first 
order phase transition [121 El 122] • 

For the numerical calculations and discussions of the nature of the phase transition we use the 
expression for the pressure at given temperature which is connected with the thermodynamical 
potential as 

P{p,T) = -[nip,T)-niO,T)], (55) 
where the thermodynamic potential has the following form 

n{p,T) = S-TS- J2 I'i^i (56) 

i=u,d,s 

and contains the internal energy S, the entropy S and the number of the ith quark Ni. 

As a matter of fact, the energy per particle of the quark system has two minimum, corresponding 
to the zeros of the pressure, the minimum at p„ = 2.25po being an absolute minimum. This implies 
that the phase transition in this model is a first order one. The results are presented in Figure 1 (left 
panel) where we have plotted the pressure in terms of Pn/po, where p„ is the neutron matter density 
and po = 0.17 fm~^ is the normal nuclear density. 

At finite temperature we conclude that for T < 56 MeV we still have a first order phase transition, 
but for T > Tci = 56 MeV the phase transition is of second order. 
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Figure 2: Masses of the u and s quarks as function of the density and temperature. 

Below this temperature there is a range of densities (p„ < 2.25po) where the pressure and/or the 
compressibihty are negative meaning that the system is in a mixing phase: a low density (p^) phase 
with massive quarks and a high density {p^) phase with light quarks of (partially) restored chiral 
symmetry (see also [T5 | IT6 | l22]). 

With this parameterization [22], — 0, and at zero temperature, the model may be interpreted 
as having a hadronic phase - droplets of light u , d quarks with a density p^*" = 2.25po surrounded by 
a non-trivial vacuum - and, above the critical density, a quark phase with partially restored chiral 
symmetry. 

The right panel in Figure 1 shows the combined effect of the temperature and density dependencies 
of the pressure. 

As the temperature increases the pressure becomes positive but the compressibility is still nega- 
tive. At T = 56 MeV and p = 1.53po, the compressibility has only on zero and we identify this point 
as the critical endpoint, the connects the first order and second order phase transition regions. 

5 77 — 77' behavior in the medium 

Figure 2 shows the temperature and density dependencies of constituent quark masses Mj (left 
panel: M„ = M^, right panel: M,). As we see at low temperature masses are almost constant and 
then they decrease with increasing temperature. The same situation can be observed in the density 
direction. 

The combined effect of density and temperature shows, as it was expected, that chiral symmetry 
is partially restored for the light quarks and is badly restored in the strange sector. There is a clear 
separation between the regions of restored and chiral symmetry in the diagram of the left panel of 
Figure 2, but the same does not happen in the right panel. 

In what concerns the mesonic behavior, before discussing our new findings, we reproduce here 
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Figure 3: a) Temperature dependence of the vr (dashed hne), rj (sohd hne) and r]' (dotted hne) 
masses. The curve 2Mu shows the temperature dependence of the quark threshold. Respective Mott 
temperatures Tmtt and Tmt] are indicated; b) Density dependence of r] (sohd hne) and rj' (dotted hne) 
masses. 

the behavior at finite temperature and zero density obtained, for instance in [23] ■ Figure 3 a) shows 
the temperature dependence of the M^r, and M^j' mesons with vanishing chemical potential. For 
comparison, the curve of 2Mu is also indicated. One see that at low temperature the M^r, masses 
are lower than masses of their constituents. In this case the integrals /2*(-Po; T) are real. The crossing 
of the TT and rj lines with 2Mu line indicates the respective Mott transition temperature Tm for these 
mesons. For our set of parameters the mesons become unbound at Tm ~ 200 MeV. One can also see 
that Tmti < Tmw Above the Mott temperature we have taken into account the imaginary parts of 
the integrals Il^iPoiT) and used a finite width approximation I23j. As we mentioned before the 
1]' - meson is unbound. 

We start the discussions of the r] and rj' in matter from the properties of dynamical quark masses 
and chemical potentials. They are plotted in Figure 4 as functions of p„. 

An important point to be noticed is that, due to the 't Hooft contribution in the gap equations, 
the mass of the strange quark decreases smoothly and becomes lower than the chemical potential 
at densities above 3.9po; which we will denote from now on as p'^. A more pronounced decrease is 
then observed, which is due to the presence of s quarks in this regime. As it will be discussed below, 
this fact is related to a change in the behavior of different observables, as compared to the region of 
lower densities. At p = p'^ (Figure 4) the difference between the chemical potentials pd = Ps has a 
maximum and, as the density increase there is a tendency to the restoration of flavor symmetry j21j . 

Now, we analyze the results for the masses of rj and rj' . 

Having in mind the constraints (j53|) and (j54p we have calculated the density dependence of 
and Mj^i masses. The results are shown on Figure 3-b. 

In order to understand these results, it is useful to study the limits of the Dirac sea of this mesons. 
They can be obtained by looking the limits of the regions of poles in the integrals /2*(Po) fj45|) and are 
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Figure 4: Constituent quark masses (solid lines) and chemical potentials, i^s — fJ-u, l-i'e (dashed lines) 
at T = MeV. 

plotted in Figure 3 (dashed point lines): Uu = 2^^, Wd = '^l^d and cUg = (At finite temperature, 
we will have a;„ = SM^, a;^ = 2Ma and ujs = 2Ms). 

As we can see, the rj' - meson lies above the quark - antiquark threshold for p„ < 2.5po and it is 
a resonant state. After that density, the rj' becomes a bound state. 

Concerning the masses of r] and rj' they exhibit a tendency to became degenerate but there are 
no clear indications of restoration of Ui{A) symmetry. In neutron matter with beta equilibrium the 
systems shows a tendency to the restoration of flavor symmetry, which is related with the presence 
of strange quarks in the medium that occurs at about p ~ 3.8po- 

Concerning the behavior in a hot and dense medium, a feature to be noticed is that the r) 
meson, that has a Goldstone boson like nature, shows more clearly the difference between the chiral 
symmetric and asymmetric phase, as shown in Figure 5. In conclusion, in the framework of our 
model, if the behavior of rj' at high temperatures and zero density suggests a possible restoration of 
the C/a(1) symmetry, the same does not happen when we explore the high density region of a dense 
or hot and dense medium. 
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Figure 5: Left panel: rj mass as function of the density and temperature; Right panel: r]' mass as 
function of the density and temperature. 



References 

[1] F. Karsch and E. Laermann, Phys. Rev. D 50 (1994) 6954; 
K. Kanaya, Prog. Theor. Phys. Sup. 129 (1997) 197. 



[2] 
[3] 
[4] 
[5] 

[6: 

[7] 

[8 
[9 



C. Roland, (PHOBOS Collaboration), Nucl. Phys. A 698 (2002) 54. 
C. Lourengo, Nucl. Phys. A 698 (2002) 13. 

Z. Fodor and S. D. Katz, Phys. Lett. B 534 (2002) 87, |hep-lat/0104001| |h^-lat/0204029] 

E. Shuryak, Comm. Nucl. Part. Phys. 21 (1994) 235. 

J. Kapusta, D. Kharzeev and L. McLerran, Phys. Rev. D 53 (1996) 5028. 

Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (1961) 345; Phys. Rev. 124 (1961) 246; 
M. K. Volkov, Ann. Phys. 157 (1984) 282. 



G. A. Sokol, L. N. Pavlyuchenko, |nucl-ex/0111020| (2001). 



K. Ohnishi, K. Fukushima, K. Ohta, Phys. Rev. C 63 (2001) 045203; Phys. Lett. B 514 (2001) 
200. 



14 



[10] U. Vogl and W. Weise, Prog. Part. Nucl. Phys. 27 (1991) 195; 
S. P. Klevansky, Rev. Mod. Phys. 64 (1992) 649; 
T. Hatsuda and T. Kunihiro, Phys. Reports 247 (1994) 221; 
P. Rehberg, S. P. Klevansky, Ann. Phys. 252 (1996) 422. 

[11] D. Blaschke et al., Nud. Phys. A 592 (1995) 561. 

[12] M. C. Ruivo, C. A. de Sousa, B. Hiller and A. H. Bhn, Nucl. Phys. A 575 (1994) 460. 
[13] M. C. Ruivo and C. A. Sousa, Phys. Lett. B 385 (1996) 39. 

[14] C. A. Sousa and M. C. Ruivo, Nucl. Phys. A 625 (1997) 713. W. Broniowski and B. Hiller, 
Nucl. Phys. A 643 (1998) 161. 

[15] M.C Ruivo, P. Costa and C. A. Sousa, in Quark Matter in Astro and Particle Physics, Eds. D. 
Blaschke, Rostock University, 2001, p. 218, hep-ph/0109234 

[16] M. C. Ruivo, C. A. Sousa and C. Providencia, Nucl. Phys. A 651 (1999) 59. 

[17] E. Witten, Nucl. Phys. B 156 (1979) 269. 

[18] G. Veneziano, Nucl. Phys. B 159 (1979) 213. 

[19] W. Broniowski and B. Hiller, Nucl. Phys. A 643 (1998) 161. 

[20] M. C. Ruivo, Hadron Physics - Effective Theories of Low Energy QCD, ed. A. Blin et al, AIP, 
New York, 2000, p. 237. 

[21] P. Costa and M. C. Ruivo, Europhysics Letters 60(3) (2002) 356. 

[22] M. Buballa and M. Oertel, Nucl. Phys. A 642 (1998) 39c; Phys. Lett. B 457 (1999) 261. 
[23] P. Rehberg, S. P. Klevansky and J. Hiifner, Phys. Rev. C 53 (1996) 410. 



15 



